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Abstract 

I set out the 2PI formalism for quantum fields in a Friedmann-Robertson- Walker Universe. I 
show how one can solve self-consistently for the quantum field evolution and the evolution of 
the cosmological scale factor, using a renormalised semi-classical Friedmann equation. In a 2PI 
coupling expansion truncated at 0{X^), I perform sample calculations demonstrating kinetic 
and chemical equilibration in this context. 



1. Introduction and setup 

Direct numerical calculations in quantum field theory often rely on lattice discretization 
of certain equations of motion, which are then solved in real time. Although this is 
straightforward in Minkowski space, where the lattice is static, when transferred to an 
expanding Universe, a number of issues need to be resolved. Firstly, the scale factor of 
the Friedmann-Robertson- Walker Universe now multiplies the lattice spacing (the lattice 
coordinates are comoving), which means that discretization errors increase in time (the 
available momentum modes redshift out of the box) . Given a finite lattice to start with, 
this imposes a cap on the amount of expansion and hence the time scales one can hope to 
reliably simulate. Secondly, the dynamics of the scale actor is determined by the energy 
momentum tensor of matter through the Friedmann equation, a quantity which in a 
quantum theory must be renormalised. 

The 2PI formalism for quantum fields in real time has been very successful recently 
in describing the approach to equilibrium and thermalisation (see for instance |1I2| V It 
is ideally suited to study strongly out of equilibrium phenomena like post-inflationary 
(p)reheating and phase transitions in the early Universe. However, unless we can argue 

Email address: ander s . tr anberg{Soulu . f i (Anders Tranberg). 



Preprint submitted to Elsevier 



16 October 2008 



that expansion is negligible (and sometimes we can), we need to extend the formalism 
to include this expansion. 

The 2PI equations of motion turn out to generalise very simply, and although an exact 
renormalisation is cumbersome, for practical numerical purposes there is a well defined 
and systematically improvable procedure by which one can renormalise at finite lattice 
spacing to the required accuracy. 

At the end of day the issue of redshift is a question of having a large enough lattice on 
a large enough computer, and can not be avoided. The approach here is to instead bypass 
the problem by only considering the phenomena, for which the expansion range fits on 
available lattice sizes. A back-of-the envelope estimate yields that a{t stop) / a{t start) < 
nx/27:, where a{t) is the scale factor, and rix is the linear scale of the lattice. This 
essentially rules out direct studies of inflation, but since T oc a~^, it allows a large range 
in temperatures and therefore most post-inflationary phenomena. We will study a single 
self-interacting scalar field in an FRW Universe with scale factor a{t), with the action 



S = J dt(fyia^{t) 



(1) 



The scale factor is treated as an external field, and evolves according to the semi-classical 
Friedmann equation SM^^H'^it) = (T°"(<)}ron, where H{t) = |[|y is the Hubble rate, and 
the energy momentum T^^ is renormalised relative to the adiabatic vacuum (see below). 
The 2PI equations of motion can then be derived in the usual way. Going to conformal 
time 77, dt = a{r])dr], and rescaling the field (/)(x, t)/ a{T) = '/5(x, <), we have 



[dl + e + M2(^)] f{t^, tI, k) = - y dT|'Y.p{^, tI', k)Fir,", 77', k) 


v' 

+ 1 dTj"j:p{r^,Tj",k)p{Tj",Tj',k), (2) 



[d^, + k^ + M\^)] p(77,,/,k) = - / VSp(77,r;",k)p(7y",77',k), (3) 



[df^+M^irj)] m = - j dri' j d^x I]^(7?, ,7', x)<^(7/), (4) 


in terms of the mean field </)(??) = {(t>{x)) and the two-point function {T(f){x)4>{y)) — 
'Piv)4>iv') = F iv 1 v' T ~ y) " ^PiViV' ~ y) ^^S^\{v ~ V')i ■ The precise form of the 
self-energy components M^, M^, E^, E^?, depends on the truncation of a diagram 
expansion, in this case a coupling expansion to C'(A^), and can be found in |f I2l3j . In 
the rescaled field, conformal time variables, the expansion of the Universe only enters 
through an effective modification of the masses + o?{'i^)w? where a" ~ 

d'^a/drf', the trivial rescaling of the momenta /cphys — k/a{ri), and the fact that time is 
stretched through t — to = a{vi)dri (see also U for related or partial implementations of 
expanding backgrounds). Renormalisation of the energy momentum tensor is ambiguous, 
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Fig. 1. Left: The effective temperature T = a;^/ ln{l/nit + 1) for each momentum mode k at different 
times, assuming a Bose-Einstein spectrum n^. Shown arc the cases of a free masslcss (black) and massive 
(red) field. The massive case is driven away from equilibrium, and an effective chemical potential is 
generated. Time increases downwards, and the physical momentum range shrinks due to redshift. Right: 
iOi^/T vs. k for an interacting massive field (black) compared to Bose-Einstein fits (red), which allow us 
to extract effective temperatures and chemical potentials. 
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Fig. 2. Effective temperatures (upper) and chemical potentials (lower) in time for an interacting massive 
field, initially in equilibrium, for different expansion rates. Blue is fastest, black slowest. Temperature 
decreases as expected, but the chemical potentials increase due to the expansion as for a free field (dashed 
lines). Only in the slowest expansion case do the interactions manage to partially compensate for the 
expansion, chemically equilibrating back towards zero chemical potential. 



and requires specification of the vacuum with respect to which we renormahse. Here we 
choose the vacuum defined by solving the vacuum 2PI-L0 (Hartree) equations as a 
function of a(t). In practice we solve them to second order in a WKB approximation 
and to second order in derivatives of a{t) (order tl). We can then define a mass 
and an energy countcrtcrm, whereby we cancel quartic and quadratic divergences, while 
neglecting logarithmic ones, at 2PI-L0. By going to higher order in WKB and H one 
can take care of the log divergences as well, but this is beyond our scope. Also, diagrams 
at NLO (A^) introduce new divergences of all orders, but these arc suppressed by the 
coupling. Details of this can be found in |5l3j . In practice, the energy momentum is 
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Fig. 3. The particle number vs. momentum for different time, starting from a way-out-of equilibrium 
state (black). Both for slow (left) and fast (right) expansion, interactions are strong enough to kinetically 
equilibrate to a final Bose- Einstein distribution (blue line). Again, we can extract effective temperatures 
and chemical potentials, to find that only at the slowest expansion rates does chemical equilibration take 
place. With fast expansion, a finite chemical potential persists. 

renormaliscd to zero in the vacuum with a relative accuracy of 10^^ for the simulation 
parameters used below. 

2. Sample applications 

We performed sample simulations using = 32'^ lattices starting at a(0)m = (mass- 
less) and a(0)m — 0.2 (massive), expanding by a factor of order 10. The coupling was 
A = 6, for which the coupling expansion is known to converge, and the LO renormali- 
sation approximation is reasonable. The expansion rate is encoded through the Planck 
mass Mpi, with larger Planck mass leading to slower expansion. Applications include: 
Free field originally in equilibrium in an expanding background (Fig. [H left). The same 
for an interacting field (Fig. [1] right and [5]). An interacting field initially far from equi- 
librium (Fig. [3]). More applications can be found in [3], including an oscillating free mean 
field (inflaton) and resonant preheating. 
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